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The analysis of the recent charged-current neutrino-nucleus scattering cross sections measured 
by the MiniBooNE Collaboration requires relativistic theoretical descriptions also accounting for 
the role of final-state interactions. In this work, we evaluate differential antineutrino-nucleus cross 
sections with the relativistic Green's function model, where the final-state interactions are described 
in the inclusive scattering consistently with the exclusive scattering using a complex optical poten- 
tial. The sensitivity to the parameterization adopted for the phenomenological optical potential is 
discussed. The predictions of the relativistic Green's function model are compared with the results 
of different descriptions of final-state interactions. 
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I. INTRODUCTION 

The MiniBooNE Collaboration has recently re- 
ported [l[ a measurement of the charged-current 
quasielastic (CCQE) flux-averaged double-differential 
muon neutrino cross section on 12 C in an energy range 
up to w 3 GeV. The neutrino-nucleus CCQE reaction 
in MiniBooNE may be considered as scattering of an in- 
cident neutrino with a single nucleoli bound in carbon, 
but it can also be sensitive to contributions from collec- 
tive nuclear effects, whose clear understanding is crucial 
for the analysis of ongoing and future neutrino oscillation 
measurements [H-Q]- 

When a dipole dependence on the four-momentum 
transferred squared Q 2 is assumed for the axial form 
factor, the nucleon axial mass Ma has been used as a 
free parameter within the relativistic Fermi gas (RFG) 
model [H [j|. Recent CCQE measurements 0, HJ re- 
ported values of Ma ~ 1-2 GeV/c 2 , significantly larger 
than the world average value from the deuterium data of 
Ma = 1-03 GcV/c 2 [TO, El- In agreement with these re- 
sults, the MiniBooNE collaboration reported values of 
M A = 1.35 ± 0.17 GcV/c 2 for the CCQE measure- 
ments Q and M A = 1.39 ±0.11 GeV/c 2 for the neutral- 
current elastic (NCE) data Q. A recent application of 
analyticity and dispersion relations to the axial vector 
form factor to find constraints for the axial mass param- 
eter using the CCQE data from MiniBooNE is presented 
in Ref. [lj and produces a value of M A = 0.85lo 22 ±0.09 
GeV, which significantly differs from the RFG extraction. 

The energy region considered in the MiniBooNE ex- 
periments, with neutrino and antincutrino energy up to 
~ 3 GeV and average energy of the muon neutrino (an- 
tineutrino) flux ps 0.79 (0.66) GeV [Hj], requires the use 
of a relativistic model, where not only relativistic kine- 
matics should be considered, but also nuclear dynamics 
and current operators should be described within a rel- 
ativistic framework. From the comparison with electron 
scattering data it is known that the RFG, although able 



of getting the basic shape and size of the response, turns 
out to be a too naive model to correctly account for im- 
portant details of the nuclear dynamics. Thus, the larger 
axial mass needed by the RFG could be considered as 
an effective value to incorporate nuclear effects into the 
calculation rather than a clear signal of a modified axial 
mass. 

At intermediate energy, quasielastic (QE) electron 
scattering calculations [1J, [l5| , which were able to suc- 
cessfully describe a wide number of experimental data, 
can provide a useful tool to study neutrino-induced pro- 
cesses. Several theoretical models have been applied in 
recent years to j/-nuclcus scattering reactions and some 
of them have been compared with the MiniBooNE data, 
both in the CCQE and in the NCE channels. At the level 
of the impulse approximation (IA), models based on the 
use of a realistic spectral function [ill [Tt}, which are 
built within a nonrelativistic framework, underestimate 
the experimental CCQE and NCE cross sections unless 
Ma is enlarged with respect to the world average value. 
The same results are obtained by models based on the 
relativistic IA (RIA) [l^-[2Cj|. However, the reaction may 
have significant contributions from effects beyond the IA 
in some kinematic regions where the experimental neu- 
trino flux has significant strength. For instance, in the 
models of Refs. (2l| - l25j the contribution of multinucleon 
excitations to CCQE scattering has been found sizable 
and able to bring the theory in agreement with the ex- 
perimental MiniBooNE cross sections without increasing 
the value of Ma- 

The role of processes involving two-body currents com- 
pared to the IA models has been discussed in Refs. fvH 
I26l428j . A careful evaluation of all nuclear effects and of 
the relevance of multinucleon emission and of some non- 
nucleonic contributions [29| - [32l | would be interesting for a 
deeper understanding of the reaction dynamic. However, 
fully relativistic microscopic calculations of two-particle- 
two-holc (2p-2h) contributions are extremely difficult and 
may be bound to model dependent assumptions. For 
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instance, the part of the 2p-2h excitations which may 
be reached through two-body meson-exchange currents 
(MEC), in particular the contribution of the vector MEC 
in the 2p-2h sector, evaluated in the model of Ref. [33| . 
has been incorporated in a phenomenological approach 
based on the superscaling behavior of electron scatter- 
ing data [H, HI [H|. The effects of MEC are impor- 
tant relative to the QE contribution, especially for the 
antineutrino cross section, where the destructive vector- 
axial interference term reduces the pure QE contribution 
and the MEC have a more significant role [34| . 

Within the QE kinematic domain, the treatment of 
the final-state interactions (FSI) between the ejected nu- 
cleon and the residual nucleus is an essential ingredi- 
ent for the comparison with data. The relevance of 
FSI has been clearly stated in the case of exclusive 
(e, e'p) processes, where the use of complex optical po- 
tentials in the distorted wave impulse approximation 
(DWIA) is required 0, QJ, • However, the pure 

DWIA approach, which is based on the use of an ab- 
sorptive complex potential, would be inconsistent in the 
analysis of inclusive scattering, where all final-state chan- 
nels should be retained and the total flux, although re- 
distributed among all possible channels due to FSI, must 
be conserved. Different approaches have been used to de- 
scribe FSI in relativistic calculations for the inclusive QE 
electron- and neutrino- nucleus scattering [42H52j ]. In the 
relativistic plane- wave impulse approximation (RPWIA), 
FSI arc simply neglected. In another approach, FSI arc 
included in DWIA calculations where the final nucleon 
state is evaluated with real potentials, either retaining 
only the real part of the relativistic energy-dependent 
complex optical potential (rROP), or using the same rel- 
ativistic mean field potential considered in describing the 
initial nucleon state (RMF). Although conserving the 
flux, the rROP is unsatisfactory from a theoretical point 
of view, since it relies on an energy-dependent potential, 
which reflects the different contribution of open inelastic 
channels for each energy, and under such conditions dis- 
persion relations dictate that the potential should have a 
nonzero imaginary term [53| . On the other hand, in the 
RMF model the same strong energy-independent real po- 
tential is used for both bound and scattering states. It 
fulfills the dispersion relations (53j and also the continu- 
ity equation. 

In a different description of FSI relativistic Green's 
function (RGF) techniques [H, |H M, M, H3 are 
used. In the RGF model, under suitable approxima- 
tions, which are basically related to the impulse approx- 
imation, the components of the hadron tensor are writ- 
ten in terms of the single particle optical model Green's 
function, whose self-energy is the Feshbach optical po- 
tential. The explicit calculation of the single particle 
Green's function can be avoided by its spectral repre- 
sentation, which is based on a biorthogonal expansion in 
terms of a non Hcrmitian optical potential % and of its 
Hermitian conjugate Calculations require matrix el- 
ements of the same type as the DWIA ones for the case of 



exclusive (e, e'p) processes in Ref. [33], but involve eigen- 
functions of both T-L and H\ where the imaginary part 
gives in one case an absorption and in the other case a 
gain of flux. This formalism allows us to reconstruct the 
flux lost into nonelastic channels in the case of the inclu- 
sive response starting from the complex optical potential 
which describes clastic nucleon-nucleus scattering data. 
Thus, it provides a consistent treatment of FSI in the 
exclusive and in the inclusive scattering and gives a good 
description of (e, e') data [H,[5l|. Due to the analyticity 
properties of the optical potential, the RGF model ful- 
fills the Coulomb sum rule [46], [53], [54| . In addition, the 
RMF and RGF reproduce also the behavior of the scaling 
function extracted from the electron scattering data [5l[ . 

These different descriptions of FSI have been compared 
in [5l| for the inclusive QE electron scattering, in Ref. 52 1 
for the CCQE neutrino scattering, and in Refs. [56|, 57 j 
with the CCQE and NCE MiniBooNE data. In Ref. [56 1 
both the RMF and the RGF give a good description of 
the shape of the CCQE experimental data and, more- 
over, the RGF can give cross sections of the same mag- 
nitude as the experimental ones without the need to in- 
crease the value of Ma- Similar results are obtained in 
Ref. [57| , where the RGF results and their interpretation 
in comparison with the NCE data from MiniBooNE are 
discussed. 

In this paper different relativistic descriptions of FSI 
for CCQE //-nucleus reactions are discussed and results 
for the double-differential cross section averaged over the 
Dfj, MiniBooNE flux are presented. The MiniBooNE col- 
laboration has accumulated an extensive data set of v ^ 
events, but it has also measured CCQE events. The 
analysis of the antineutrino data is currently ongoing (58j 
and some preliminary results can be found in the Mini- 
BooNE website [5{|. When available, the antineutrino 
measurements will be an additional source of informa- 
tion about the weak charged-current lepton-nucleus in- 
teraction and, combined with the corresponding neutrino 
data, will provide important insight about the role of the 
longitudinal, transverse, and interference vector-axial re- 
sponses which enter the cross sections. Indeed, being 
aware of the interpretative questions which may be con- 
nected to the fact that the neutrino and the antineutrino 
fluxes at MiniBooNE are different, with the flux signif- 
icantly smaller and with lower average energy than the 
one, measurements of both reactions would be use- 
ful to clarify the role of nuclear effects in the analysis of 
lcpton-nuclcus processes. 



II. RESULTS AND DISCUSSION 

In all the calculations presented in this work the 
bound nucleon states are taken as self-consistent Dirac- 
Hartree solutions derived within a relativistic mean field 
approach u sing a Lagrangian containing a, uj, and p 
mesons [60 62]. In the RGF calculations we have used 
two parameterizations for the relativistic optical poten- 
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Figure 1. (color online) Flux-averaged double-differential 
cross section per target nucleon for the CCQE 12 C(^ M ,/i + ) 
reaction as a function of T M for four angular bins of cosi?^ 
calculated with the RGF-EDAD1 (solid lines) and the RGF- 
EDAI (dashed lines). The dotted lines are the rROP results 
calculated with the EDAI potential and the dot-dashed lines 
are the RPWIA results. 



tial: the Energy-Dependent but A-Indcpcndcnt EDAI 
(where the A represents the atomic number) and the 
Energy-Dependent and A-Dependent EDAD1 complex 
phenomenological potentials of Refs. [63T465T ] . which are 
fitted to elastic proton scattering data in an energy range 
up to 1040 MeV. We note that whereas EDAD1 is a 
global parameterization, EDAI is a single-nucleus param- 
eterization, which is constructed to fit scattering data 
just on 12 C and, as such, does have an edge in terms 
of better re pro duction of the elastic proton- 12 C phe- 
nomenology [63[ compared to EDAD1, and also leads to 
a better description of the inclusive quasielastic (e, e') 
cross sections, as well as to CCQE and NCE results in 
better agreement with the MiniBooNE data [5(1 [57| ■ 

In Fig. [T] the CCQE double-differential 12 C(^,/i+) 
cross section per target nucleon integrated over the Mini- 
BooNE flux is shown as a function of the muon kinetic 
energy T M for four angular bins of cos i9 M , where fi^ is the 
muon scattering angle, ranging from forward to backward 
angles. In the RPWIA calculations FSI are completely 
neglected. The rROP results, where calculations are per- 
formed with a pure real optical potential, are usually 15% 
lower than the RPWIA ones. The rROP generally under- 
estimates the v experimental data unless a larger axial 
mass, e.g., Ma ~ 1.3 — 1.4 GeV/c 2 , is used. However, in- 
dependently of its comparison with the data, the rROP 
model, which is based on an energy-dependent poten- 
tial, has important physical drawbacks [46l l52l . l53l . [56j . 
The RGF cross sections with both optical potentials are 



larger than the RPWIA and the rROP ones. The differ- 
ences between the RGF results with the two optical po- 
tentials are clearly visible. For instance, the EDAI and 
EDAD1 potentials yield differences by about 15% — 20% 
to the cross section in the peak region for the forward an- 
gle scattering bins cos# M = 0.85, 0.75. Somewhat closer 
predictions are obtained in the bin cosfJ^ = 0.45, while 
for the backward angular bin cost?^ = —0.15 the dif- 
ferences are enhanced up to 25%, but the magnitude of 
the cross sections is significantly reduced. We note that 
the relative differences between the RGF results with the 
two optical potentials are somewhat larger in neutrino 
scattering j56[. The different behavior of the RGF in 
neutrino and antincutrino scattering is related to the rel- 
ative strength of the vector-axial response, which is con- 
structive in v scattering and destructive in v scattering 
with respect to the longitudinal and transverse ones [52j . 
Moreover, the differences between the neutrino and an- 
tineutrino MiniBooNE fluxes, make the comparison be- 
tween the results of ^-nucleus and P M -nucleus scattering 
not straightforward. 

The comparison between the RGF results obtained 
with the EDAI and EDAD1 potentials can give an idea 
of how the predictions of the model are affected by un- 
certainties in the determination of the phenomenological 
optical potential. The differences depend on the energy 
and momentum transfer and arc essentially due to the 
different values of the imaginary parts of the two po- 
tentials, which account for the overall effects of inelas- 
tic channels and are not univocally determined from the 
elastic phenomenology. In contrast, the real terms are 
very similar for different parameterizations and give very 
similar results. In the rROP calculation shown in Fig. [TJ 
the real part of the EDAI potential has been used, but a 
calculation with EDAD1 would give in practice the same 
result. The results in Fig. [T] stress the importance of FSI 
and, in particular, of the imaginary part of the relativis- 
tic optical potential, which plays a different role in the 
different approaches. In the rROP, the imaginary part is 
neglected and the total flux is automatically conserved. 
The RGF results presented here contain the contribu- 
tion of both terms of the hadron tensor in Eq. (61) of 
Ref. (46[. The calculation of the second term, which is 
entirely due to the imaginary part of the optical poten- 
tial, is a hard and time consuming numerical task which 
requires the integration over all the eigenfunctions of the 
continuum spectrum of the optical potential. Numerical 
uncertainties on this term are anyhow under control and, 
from many calculations in different kinematics, have been 
estimated at most within 10%. In the RGF, the imagi- 
nary part redistributes the flux in all the final-state chan- 
nels and, in each channel, the loss of flux towards other 
inelastic channels (either multinucleon emission or non- 
nucleonic excitations) is compensated for the inclusive 
scattering making use of the dispersion relations. The 
larger cross sections in the RGF arise from the trans- 
lation to the inclusive strength of the overall effects of 
inelastic channels. These contributions are not included 
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Figure 2. (color online) Flux-averaged double-differential 
cross section per target nucleon for the CCQE 12 C(^ M ,/i + ) 
reaction as a function of cosi9 M for two bins of T M calcu- 
lated with the RGF-EDAD1 (solid lines) and the RGF-EDAI 
(dashed lines). The dotted lines are the rROP results calcu- 
lated with the EDAI potential and the dot-dashed lines are 
the RPWIA results. 



Figure 3. (color online) Total CCQE cross section per target 
nucleon as a function of the neutrino energy E v (upper panel) 
and of the antineutrino energy Ep (lower panel) calculated 
with the RGF-EDAD1 (solid lines), the RGF-EDAI (dashed 
lines), the rROP (dotted lines), and the RPWIA (dot-dashed 
lines). The experimental data for neutrino scattering are from 
MiniBooNE 



explicitly in the RGF model, but can be recovered, at 
least to some extent, in the RGF by the imaginary part 
of the phcnomcnological optical potential which reincor- 
porates these processes in the reaction. 

In Fig. [2] the CCQE flux-averaged double-differential 
12 C(P M ,/x + ) cross section per target nucleon is displayed 
as a function of cos i? M for two bins of T M . Similar consid- 
erations can be made for these results: the rROP gives 
the lower results and the RGF cross sections with both 
optical potentials are larger than the RPWIA and the 
rROP ones. All the models, however, tend to produce 
similar results for backward scattering angles, where the 
cross sections are sensibly reduced. The differences be- 
tween the RGF results with the two optical potentials 
are visible but they are somewhat reduced with respect 
to the correspon ding calculations for neutrino scattering 
in Fig. 2 of Rcf. 56[ . This is essentially due to the com- 



bined effects of the differences between the and the 
fluxes and of the destructive contribution of the vector- 
axial interference response. 

In Fig. |3] the total QE cross sections per target nucleon 
for neutrino and antineutrino scattering are displayed as 
a function of the neutrino (antineutrino) energy E v (Ep). 
The neutrino results, which have been already presented 
in Ref. jfjfjj . are compared with the experimental data 
from MiniBooNE In Refs. [H H3 it is shown that 
DWIA models where FSI effects arc accounted for by 
means of real potentials, like rROP and RMF, produce 
similar results which all underestimate the total CCQE 



MiniBooNE cross section, whereas the inclusion of 2p-2h 
MEC enhances the results. Larger cross sections are ob- 
tained in the RGF with both optical potentials. The 
differences between RGF-EDAI and RGF-EDAD1 arc 
clearly visible, RGF-EDAI being in good agreement with 
the shape and magnitude of the experimental cross sec- 
tion for neutrino scattering. We observe that EDAI is a 
single-nucleus parameterization, which is constructed to 
better reproduce the elastic proton- 12 C phenomenology, 
and gives electron scattering cross sections in fair agree- 
ment with the experimental data. We also note that the 
antineutrino cross section does not saturate in the energy 
range up to 2.5 GeV which we have considered. Also 
in this case, the differences between EDAI and EDAD1 
are due to the different imaginary parts of the two po- 
tentials, particularly for the energies considered in kine- 
matics with the largest T M . These kinematics give large 
contributions to the total cross section and enhance the 
differences between the two RGF results. 



III. SUMMARY AND CONCLUSIONS 

In this paper, we have compared the predictions 
of different relativistic descriptions of FSI for CCQE 
antineutrino-nucleus scattering in the MiniBooNE kine- 
matics. In the RPWIA, FSI arc simply neglected; in the 
rROP, they arc described retaining only the real part of 
the relativistic energy-dependent optical potential; in the 
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RGF, the full complex optical potential, with its real and 
imaginary parts, is used to account for FSI. All final-state 
channels are included in the RGF, where the flux lost in 
each channel is recovered in the other channels by the 
imaginary part of the optical potential making use of the 
dispersion relations and the total flux is conserved. The 
RGF is able to give a good description of the (e, e') data 
in the QE region and it is also able to describe both the 
MiniBooNE CCQE and NCE neutrino data without the 
need to change the value of the axial mass. 

The enhancement of the RGF cross sections compared 
to the cross sections calculated with other descriptions 
of FSI is obtained also in the case of antineutrino scat- 
tering. The larger results of the RGF can be ascribed 
to the contribution of reaction channels which are not 
included in the other models. For instance, rescattering 
processes of the nucleon in its way out of the nucleus, 
non-nucleonic A excitations, which may arise during nu- 
cleon propagation, with or without real pion production, 
or also multinucleon processes. These contributions are 
not included explicitly in the RGF model, but can be 
recovered, to some extent, by the imaginary part of the 
phcnomcnological optical potential. We cannot disentan- 
gle the role of different reaction processes and explain in 
detail the origin of the recovered strength. It would be 
anyhow interesting for a comparison to disentangle in the 
phenomenological optical potential the contributions due 
to non-nucleonic inelasticities and extract a "purely nu- 
clconic" optical potential that could be used in the RGF 
approach. 

The RGF predictions are also affected by uncertain- 
ties in the determination of the phcnomcnological opti- 
cal potential, which is not univocally determined from 
the elastic phenomenology. In the case of antineutrino 
scattering the differences between the RGF cross sections 
calculated with the EDAI and EDAD1 optical potentials 
are somewhat reduced with respect to the correspond- 
ing cross sections calculated for neutrino scattering. The 
differences between the RGF results for neutrino and an- 
tineutrino scattering may be ascribed to the different 
and fluxes and to the strength of the vector-axial re- 
sponse. A better determination of the phenomenological 
relativistic optical potential, which closely fulfills the dis- 



persion relations, would reduce the theoretical uncertain- 
ties on the RGF results and, therefore, deserves further 
investigation. 

The analysis of MiniBooNE CCQE and NCE data with 
theoretical models based on the IA and including only 
one-nucleon knockout contributions usually requires a 
larger value of Ma to reproduce the magnitude of the ex- 
perimental cross sections. The calculations required for 
the theoretical analysis must consider the entire range 
of the relevant MiniBooNE neutrino energies and ad- 
ditional complications may arise from the flux-average 
procedure to evaluate the cross sections, which implies a 
convolution of the double-differential cross section over 
the neutrino spectrum. Because of uncertainties associ- 
ated with the flux-average procedure, the MiniBooNE 
cross sections can include contributions from different 
kinematic regions, where other reaction mechanisms than 
one-nucleon knockout are known to be dominant [l6l.ll7j|. 

Models including other contributions than one-nucleon 
knockout, like our RGF, but also the model of Refs. plT - 
[23| , where multinucleon components are explicitly in- 
cluded, are able to describe the MiniBooNE neutrino 
data without the need to change the value of the axial 
mass. Despite their differences, the two models seem to 
go in the same direction. In the RGF, however, the en- 
hancement of the cross section cannot be attributed only 
to multinucleon processes, since we cannot disentangle 
the role of the various contributions included in the phe- 
nomenological optical potential. In order to clarify this 
point a careful evaluation of all nuclear effects and of 
the relevance of multinucleon emission and of some non- 
nucleonic contributions would therefore be highly desir- 
able. 
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